Force induced triple point for interacting polymers 
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We show the existence of a force induced triple point in an interacting polymer problem that allows 
two zero- force thermal phase transitions. The phase diagrams for three different models of mutually 
attracting but self avoiding polymers are presented. One of these models has an intermediate phase 
and it shows a triple point but not the others. A general phase diagram with multicritical points in 
an extended parameter space is also discussed. 
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Soon after the discovery of the double hehcal struc- 
ture of double stranded DNA (dsDNA), its melting by 
change of temperature or the pH of the solvent was 
recognizedijy . Only recently it has been realized that 
there can be a force-induced unzipping transition^ of 
a DNA with a force applied solely at one end. In both 
the thermal and the forced cases, the phase transition 
takes a double stranded form to two single strands. The 
prediction of the unzipping transition based on interact- 
ing Gaussian chainS|^was immediately reconfirmed by 
a dynamical approach^. Results on the unzipping tran- 
sition are now available from extensive exact solutions 
of lattice models £], simple models of quenched- 
averaged DNA , Monte Carlo simulations of a three di- 
mensional model with self- and mutually-avoiding walks 
etc., on the theoretical front, and for real DNA from 
experiments 0. From the theoretical results, it emerges 
that the qualitative features of the unzipping transition 
are insensitive to the dimensionality (d) of the models 
and are seen even in two dimensional models. These re- 
sults include a re-entrance0,0,|H in the low temperature 
region. 

In the systems studied so far, there is only one zero 
force thermal phase transition. In such cases, with two 
intensive variables, temperature T and force g, there is 
an unzipping transition line, g — gdT), in the g ~ T 
plane, demarcating the bound or zipped phase from the 
unzipped phase. This gives the unzipping phase bound- 
ary. In case there are more than one transition, the phase 
diagram will be influenced by the intermediate phases. 
Our aim is to determine such a global phase diagram for 
the unzipping transition where the existence of an inter- 
mediate phase leads to a triple point in the g — T plane. 

We introduce three models. A, B and C, which dif- 
fer in the nature of the mutual interaction of the two 
strands. These models are general enough to be de- 
fined in any d though, for computational limitations, 
we consider the two dimensional square lattice version 
only. Let us consider two linear polymer chains which 
are mutually-attracting-self-avoiding .IQJ in nature. On 
a square lattice, the polymers are not allowed to cross 
each other as shown in Fig. ^ Monomers are the sites 
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FIG. 1: Figures (a-d) represent the possible conformations 
of models A and B. In model A any monomer of one strand 
can interact with any monomer of the other strand, (a) and 
(b) are two possible states with (c) representing a possible 
ground state. For model B, the ith monomer of one strand 
interacts with the ith monomer of the other strand. In this 
case (a) represents the ground state, (d) represents a partial 
bound state. Note that (d) differs from (b) in the nature of 
interactions represented by the dotted lines. In model B, (c) 
has no valid interaction and would represent an open state. 
Figures (e,f,g) represent model C where bases are on the links 
of the strands with short stubs representing orientation, (e) 
the completely zipped state, (f) a non pairing configuration, 
and (g) a partial bound configuration. In all cases dotted lines 
show the attractive interactions. 



occupied by the polymers and the interactions are among 
the "bases" (defined next) so that an interacting pair 
may also be called a base pair. In models A and B, the 
bases are the monomers while in C the bases are associ- 
ated with the links. We allow an attractive interaction 
between monomers or bases only if they are of opposite 
strands and are nearest neighbours on the lattice. The 
nearest neighbour interaction mimics the short range na- 
ture of the hydrogen bonds. In model A, any monomer 
of one strand can interact with any monomer of the other 
strand. In model B, monomer i of one strand can inter- 
act only with the zth monomer of the other strand. This 
model is similar to the models of DNA studied earlier j^. 
In model C, we associate the bases with the bonds of 
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FIG. 2; Schematic figures for a force (a) applied at one end, 
and (b) applied at the middle. Fixed end points are indicated 
by stars. 



the strands. In one strand the bases point towards the 
right while on the other they are on the left, as one tra- 
verses the chains seauentiallv|l2j. Pairing between com- 
plimentary strands takes place only when the bases ap- 
proach each other directly without the strands coming 
in-between, as shown in Fig. Such a model has been 
studied recently in the context of relative stabilities of 
DNA hairpin structures Ts^. The figure also shows a sit- 
uation where pairing cannot take place and such con- 
figurations are counted in the unbound states. The in- 
teractions are shown by dotted lines in Fig. The 
interaction energy in all the cases is taken as — e(e > 0) 
and we shall choose e — 1. In all cases the monomers at 
one end (index=l) of each strand are always kept fixed 
occupying nearest neighbour sites on the lattice. 

For model A the ground state is a spiral of the type 
shown in Fig. ^ while for model B and C it is a zipped 
state as in Fig. ^,e. Because of the constraint of hold- 
ing the monomers with index 1 at nearest neighbour 
sites, model A is also equivalent to a diblock copolymer 
modelfl^ which has two phase transitions in zero force 
with increasing temperature. The polymers go from a 
compact spiral-like phase to a zipped phase which then 
melts at a still higher temperature. It is this intermedi- 
ate phase that is of interest to us. In view of the simple 
ground state in models B and C, no such intermediate 
phases are expected or known. 

A force is apphed at one end (Figl^a)) or at the mid- 
dle (Fig. Olb)) of the chains, in the y-direction. The 
contribution to energy by this force is —gy where y is 
the absolute distance in the y-direction between the two 
strands at the point of application of the force. A re- 
cent study showed a rich phase diagram when a force 
is applied somewhere in the interior. Furthermore, such 
situations occur in many processes like gene-expression 
where RNA is formed in bubbles or eye-type configura- 
tions on DNA. Motivated by these, we consider the case 
of a force applied at the mid point. A low temperature 
analysis shows that the phase diagram is expected to be 
different from the end case. We compare the g — T phase 
diagrams of model A vs models B and C. 

The thermodynamic properties associated with the un- 
zipping transition are obtained from the partition func- 
tion which can be written as a sum over all possible con- 




FIG. 3: Variations in fluctuation of m with temperature (T) 
for g = 0.0 and 0.1 are shown for model A (a) and model B (b). 
In model B and C we do not find two peaks. The first peak 
in (a) corresponds to a transition where spiral conformation 
goes over to a zipped conformation. The second peak in (a) 
and the peak in (b) correspond to the unzipping transition. 

figurations 

Zn{u;,u) = Y,C{m,y)u:^uy . (1) 

Here TV is the chain length of each of the two strands, 
LJ = exp(l/r) is the Boltzmann weight associated with 
each base pair (taking the Boltzmann constant k-Q = 1) 
and m is the total number of intact base pairs in the 
chain. Finally u is the Boltzmann weight, exp((7/T) asso- 
ciated with force. C{m, y) is the number of distinct con- 
figurations having m base pairs whose end (or mid) points 
are at a distance y apart. We have obtained C{m,y) for 

< 16 monomers in two dimension (d = 2) and analyzed 
the partition function through series analysis. We prefer 
this technique because in this case the scaling corrections 
are correctly taken into account by suitable extrapolation 
technique. To achieve the same accuracy by the Monte- 
Carlo method, a chain of about two orders of magnitude 
larger than in the exact enumeration method has to be 
considered 0|. 

The reduced free energy per base pair is found from 
the relation 

Giiu,u) = lim ^ log ZiLO,u). (2) 

The limit A'^ — )■ oo is found by using the ratio method 
|l6l | for extrapolation. The transition point for zero 
force (i.e. thermal melting) can be obtained from the 
plot of G(w, u) versus w or from the peak value of 
d^G/d{\nu!)^ . For self-avoiding walk at m = 1, we find 
Tc = l.l±0.1,0.61±0.08,0.59±0.08for models A, B and 
C respectively. With a force (u ^ 1), the phase boundary 
is obtained from the fluctuation in m. Fig. |31 shows the 
variation of fluctuation of m with temperature for model 
A and B. 

The phase diagrams are shown in Fig. 2| 
The phase diagram for Model A has several distinctive 
features. We see two peaks in the temperature depen- 
dence of the fluctuation in m for small g. At low T there 
is a spiral structure with vanishing entropy but as T in- 
creases (T > 0.5 ) the chains acquire a zipped state. 
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FIG. 4; Force- Temperature (g-T) phase diagrams for the end 
and the mid case, (a) for model A, (b) model B and model C 
within errorbar. 



FIG. 5: Possible phase diagrams without any triple point in 
the g — T plane. Phase 2 (a) may persist till T — under force, 
(b) may persist under force but disappears at T = 0, or (c) 
exists only over a limited range of force and temperature. In 
(a) the sequence of 1-2- unzipped transitions will be observable 
in the groundstate itself, while in (c) a reentrant behaviour 
will be observable within phase 1 itself. 



eventually unzipping at higher temperatures. The low 
temperature peak is the transition where the spiral state 
goes over to the zipped state while the second peak is 
the unzipping transition. Although model C is a more 
realistic model for representation of dsDNA, we find the 
phase diagram to be almost the same as that of model B. 
This may be taken as an la posteriori justification of us- 
ing model B for DNA. (However, in case of hetrogeneous 
chains, model C may give intermediate states by forming 
the hairpin kind of structure as shown in Figl e. Such 
structures are not possible in model B.) 

Our results are based on the TV = 16 (32 monomers) 
enumerations. For this length there is significant surface 
contribution. The ground state energy is Eq = — (2iV — 
0(\/lV)) while the ground state energy is Eq = —Ne 
for both models B and C for a DNA with N monomers 
or bases each. The 0{\/N) correction for model A 
comes from the monomers on the boundary. If we ig- 
nore the surface contribution (valid for large N) the 
spiral to zipped state transition temperature Td may 
be estimated from a simple energy balance argument. 
The spiral state has an energy Eo ~ —2N with negli- 
gible entropy while the zipped state has the free energy 
— ~N — NThifih, where In/i^ is the entropy per base 
pair of the zipped phase. Equating these two we get 

T,i = l/ln^b-0(7V-i/2), (3) 

where the surface correction has also been shown. If we 
use fib = 2.682, the connectivity constant for the square 
lattice self-avoiding walk[loj . we find T^i ~ 1.04 which 
is close to the value known from other estimates 
Since the surface correction lowers the estimate of Td, 
our value of Td = 0.5 is consistent with the form of Eq. 

m 

The surface contribution (0{\/N) in the spiral state 
for small chains is not negligible, and it helps in the sta- 
bilization of of the spiral structure over the zipped phase 
for small g because of the extra energy ~ gVN- As a 
result, the phase boundary obtained has a finite slope. 
For large A^, the force term should not affect the spiral 
to zipped transition since the force affects the relative 
separation of the strands. Therefore one would expect 



a phase boundary parallel to the g-axis meeting the un- 
zipping phase boundary at a triple point. Other possible 
phase diagrams are shown in Fig. From thermody- 
namic stability analysis 17], it is known that the angle 
between two coexistence lines at a triple point in a phase 
diagram has to be less than tt. Therefore, a discontinu- 
ity in the slope in the unzipping phase boundary in the 
g — T plane is expected at the potential triple point. The 
inset in Fig. 0^ shows the meeting of the two unzipping 
boundaries. Our main result is the prediction of the coex- 
istence of all the three phases in the g — T phase diagram 
for model A - a phenomenon not possible without force. 

At T = 0, the critical force can be found from a match- 
ing of the ground state energy with the energy of the 
completely unzipped state. If the force stretches the 
strands completely, the unzipped state has the energy 
— —2Ng taking the bond length to be unity. Com- 
paring this with the bound state energy Eq, we see 
gc{T = 0) = Ea/2N. For large N, Eq ^ 2N and so we 
get gc{T = 0) = 1, while it is 0.5 for models B and C. For 
finite N, if the surface correction is taken into account, 
then gc^l - {1/V2N) ~ 0.8232 for N = 16. The value 
shown in Fig. 2|is very close to this estimate rather than 
the large N value. Any extrapolation therefore would 
have to take into account this surface correction prop- 
erly. In absence of any large N data we do not attempt 
such extrapolations here for model A. 

The low temperature phase boundary for the mid case 
can be obtained by an extension of the T = argument 
given above. Following Ref. 0, let us consider the situa- 
tion with the force applied at a position sN, (0 < s < 1), 
from the fixed end. The unzipped state has the energy 
= —2sNg taking the bond length to be unity. Comparing 
this with the bound state energy we see gc{T 0) = a/2s 
where a = 2 for model A but a = 1 for models B and C. 
We see a factor of 2 difference in the end vs mid case. 

Close to T = 0, the free energy of a bubble of 2m base 
pairs with respect to the completely bound state is 

AF = 2m(a + Tin fib)- gm, (4) 

where In is the entropy per base pair of the bound state 
which is lost on unzipping m pairs. For the mid case a 
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FIG. 6: Proposed phase diagram in g — T — 5 space. M is 
the multicritical point (at zero force) where the zipped phase 
vanishes as 5 is made negative (attractive three-body interac- 
tion). The dotted Une in the g = plane is a crossover line 
for a 5-interaction dominated phase to go over a non-spiral 
compact phase dominated by interactions at corners. The 
triangular region in the g — T plane for the zipped phase van- 
ishes at M. M is therefore the end point of the line of triple 
points (thin line through the triangles). 



large m is favourable (m < N/2) if g > gc{T) = 2{a + 
Tin fib). In case In /it, > 0, we see a positive slope and 
hence a re-entrance, in the low temperature region. In all 
the three models studied here, In /if, ~ 0, and therefore 
no re-entrance. 

The results for model A can be extended to a case 
with a three body interaction that stabilizes the spiral 
phase[l3- Let there be an interaction 5 that favours a 
configuration with a monomer of one strand sandwiched 
between the monomers of the other strand on the two 
nearest neighbors on the same axis. Let us call such 
contact as 5-contact. For 5 < (i.e. attractive three 
body interaction) the compact spiral phase is the ground 
state as at (5 = 0. For large negative S, there is only one 
denaturation transition. Therefore a multicritical point 
occurs in the g = 0{S — T) plane. This is indicated by 
M in Fig. |3 For repulsion, 5 > 0, one may have a 
compact ground state without any 5-contact. In this case 
the ground state need not be a spiral and will have a 
mixture of double contacts (e.g. at corners) and single 
contacts. Let these occur in proportion^ : (l—p) {p < 1). 
The ground state energy is then — iV(l +p). A meanfield 
estimate of equal proportion gives E = —{3/2)N. This 
energy needs to be compared with the energy of the spiral 
state, namely -2N + 2N6 for iV oo. At T = 0, a 
change in the ground state structure occurs at (5 = 1 — 
p « 0.5. In the 5 — T phase diagram one will have a 
crossover line (dotted line in Fig. O which should end 
in the 5 — T plane at the multicritical point which is the 
confluence of the two zero-force thermal transitions. For 
large negative i5, the compact phase goes over directly 



to the unzipped phase and therefore no triple point can 
exist on the unzipping boundary. The g — T — 5 phase 
diagram is schematically shown in Fig. |H1 

To summarize, we have shown in this paper the exis- 
tence of a force induced triple point in a problem of inter- 
acting polymers where the interactions allow an entropy- 
stabilized phase before melting. This opens up the pos- 
sibility of much richer phase transitions and other multi- 
critical behaviours in polymeric systems when subjected 
to force. We hope that single molecule experiments would 
be able to explore experimentally this new area in poly- 
mers. 
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